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Abstract 



The calculation of the moments in inclusive B meson semileptonic decays is up- 
<^ graded to 0{a^). The first three moments of the lepton energy and invariant hadronic 

^^ mass distributions are computed for arbitrary cuts on the lepton energy and in var- 

*K^ ious renormalization schemes, finding in general small deviations from the O(ag/3o) 

1^ calculation. I also review the relation between MS and kinetic heavy quark masses. 
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1 Introduction 

The Vcb element of the Cabibbo-Kobayashi-Maskawa quark mixing matrix is naturally deter- 
mined in semileptonic B decays to charmed hadrons. In the case of exclusive final states, like 
for B — > D^*Hu, the corresponding form factors have to be computed by non-perturbative 
methods, for instance on the lattice. On the other hand, in the case of inclusive semileptonic 
decays B — >■ Xcip the existence of an Operator Product Expansion (OPE) ensures that 
non-perturbative effects are suppressed by powers of the bottom mass rub and are param- 
eterized by a limited number of matrix elements of local operators which can be extracted 
from experimental data. The total inclusive width and the first few moments of the kine- 
matic distributions are therefore expected to be well approximated by a double series in as 
and Aqcd/"'^;) [D E]. The general strategy for the inclusive determination of \Vcb\ consists 
in extracting the most important non-perturbative parameters, including the heavy quark 
masses, from the moments and to employ them in the OPE expression for the total width. 
A determination of \Vcb\ follows from the comparison with the experimental total rate. 

The main ingredients for an accurate analysis of the experimental data have been known 
for some time. Two implementations are currently employed by the Heavy Flavour Averaging 
Group (HFAG) [3j, based on either the kinetic scheme [H |5], El [7] (see also [8] for an earlier 
fit) or the IS* scheme [9]. They both include terms through 0(a^/3o) [lO] and 0(l/m^) [H] 
but they use different perturbative schemes and approximations, include a slightly different 
choice of experimental data, and estimate the theoretical uncertainty in two distinct ways. 
According to the latest global fits the two methods yield very close results for \Vcb\ [12] . 

The reliability of the inclusive method rests on our ability to control the higher order 
contributions in the double series and to constrain quark-hadron duality violation, i. e. effects 
beyond the OPE. The calculation of higher order effects allows us to verify the convergence 
of the double series and to reduce and properly estimate the residual theoretical uncertainty. 
Duality violation effects [13] can be constrained a posteriori, by looking at whether the 
OPE predictions fit the experimental data. This in turn essentially depends on precise 
measurements and precise OPE predictions. As the experimental accuracy reached at the B 
factories is already better than the theoretical accuracy for all the measured moments, any 
effort to improve the latter is strongly motivated. 

There has been recent progress in this direction. First, the complete two-loop perturba- 
tive corrections to the width and to the moments of the lepton energy and hadronic mass 
distributions have been computed in Refs. [HI |T5l |16]. This represents an important im- 
provement with respect to the 0(Q;g/3o) or BLM corrections, which in B decays generally 
dominate the O(a^) effects when a^ is normalized at nib- The main goal of this paper is to 
incorporate the new non-BLM corrections in the calculation of the semileptonic moments 
and to discuss their numerical impact in different schemes. 

Higher order power corrections have also been recently considered: a first analysis of 
0{l/m'^) and 0{l/mQ) effects has been presented in [17]. In the higher orders of the OPE 
there is a proliferation of operators and therefore of non-perturbative parameters: as many 
as nine new expectation values appear at 0{l/ml). Since they cannot be fitted from experi- 



ment, the authors of Ref. [TTj estimated them in the ground state saturation approximation 
and found a relatively small +0.4% effect on \Vcb\- While this sets the scale of higher order 
power corrections, it is for the moment unclear how much the result depends on the use of 
that approximation. 

Another important source of theoretical uncertainty are the ©(agAq^j-i/m^) corrections 
to the width and to the moments. Only the 0{asfi'l/'ml) terms are known [TH] at present. 
A complete calculation of these effects has been recently performed in the case of inclusive 
radiative decays [19], where the 0{as) correction increase the coefficient of /i^ in the rate 
by almost 20%. The extension of this calculation to the semileptonic case is in progress. 

Inclusive semileptonic B decays are not only useful for the extraction of Vcb'- the moments 
are sensitive to the values of the heavy quark masses and in particular to a linear combination 
of rric and nib [20], which to good approximation is the one needed for the extraction of \Vcb\ 
[12] . The b quark mass and the OPE expectation values obtained from the moments are 
crucial inputs in the determination of \Vub\ from inclusive semileptonic decays, see e.g. [21] 
and refs. therein. The heavy quark masses and the OPE parameters are also relevant for 
a precise calculation of other inclusive decays like B — )■ Xg'y [22] ■ On the other hand, ttIc 
and rub can now be measured precisely from data on charm and bottom production in e~^e~ 
annihilation [23ll2^. and from the moments of heavy quark current correlators computed on 
the lattice [25]. After checking the consistency of the constraints on rric^b from semileptonic 
moments with these precise determinations, we should therefore include them as external 
inputs in the semileptonic fits (see [12] for a first attempt), eventually improving the accuracy 
of the \Vub\ and \Vcb\ determinations. 

In this paper we set up the tools for such improved analyses, extending the code developed 
in [6] and employed by HFAG to allow for NNLO fits in arbitrary mass schemes. In particular, 
in order to avoid the non-negligible uncertainty in the conversion from the kinetic scheme 
to MS masses, we directly employ an MS definition for the charm quark mass. In addition, 
the study of scheme and scale dependence gives us useful elements for an evaluation of the 
residual theoretical uncertainty. 

It is worth recalling that in the current semileptonic fits the sensitivity to rub and /i^ 
is enhanced by the inclusion of the first two moments of the photon energy distribution 
in i? — )■ Xs7, which are now measured with good precision. Their impact on the fits is 
equivalent to that of a loose constraint on rUb with 80MeV uncertainty. However, the lower 
cut on the photon energy introduces a sensitivity to the Fermi motion of the 6-quark inside 
the B meson [7j and there are poorly known subdominant contributions not described by 
the OPE even in the absence of a photon energy cut [26]. All this makes radiative moments 
different from semileptonic ones: in the following we will concentrate on the latter. 

The paper is organized as follows: in the next Section we discuss the perturbative cor- 
rections of the lepton energy moments and the implementation of the NNLO contributions 
and we give tables of results in a few typical cases. In Section 3 we do the same for the 
moments of the invariant hadronic mass distribution. In Section 4 we review the relation 
between kinetic and MS heavy quark definitions and provide numerical conversion formulas 
with uncertainty. Section 5 briefly summarizes the main results of the paper, while the 



Appendix updates the prediction of the total rate of 5 — )■ Xf.ii' and gives an approximate 
formula for the extraction of \Vcb\- 

2 Lepton energy moments at NNLO 

In this section we consider the first few moments of the charged lepton energy spectrum in 
inclusive b — )■ ciu decays. They are experimentally measured with high precision (better 
than 0.2% in the case of the first moment) but at the S-factories a lower cut on the lepton 
energy, E^ > Ecut, is applied to suppress the background. In fact, experiments measure the 
moments at different values of Ecut- The relevant quantities are therefore 
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which is needed to relate the actual measurement of the rate with a cut to the total rate, 
from which one conventionally extracts Vd,- Since the physical information that can be 
extracted from the first three linear moments is highly correlated, it is convenient to study 
the central moments, namely the variance and asymmetry of the lepton energy distribution. 
In the following we will consider only R* and 

il = {Ei)Et>Ecut, ^2,3 = {{Ei - {El)) ' )Ei>Ecut ■ (3) 

These four observables are all functions of rrib and of the two dimensionless quantities 

r=—, e = —^ ■ 4 

nib rUb 

Our calculation follows closely the one described in [6]. Apart from the inclusion of the 
complete O(a^) corrections there are a few changes and improvements worth mentioning: i) 
the complete 0{as) and 0(q;^/3o) corrections to the charged leptonic spectrum have been first 
calculated in [27] and [28] . While they can be computed numerically for any value of ^ and 
r, a numerical integration would slow down the fitting routines significantly, hence the need 
for interpolation formulas that must be accurate in a wide range of ^ and r values, keeping 
in mind that the values of mc^b may differ considerably from one scheme to another. The 
accuracy of the interpolation formulas is important because of the cancellations that will 
be discussed shortly and has been improved wrt [6] using high precision numerical results 
based on y^. The approximations used for the 0{as) and 0(a^/3o) corrections to the linear 
moments are now quite precise and their range extends to < ^ < 0.75 and 0.18 < r < 0.28; 



ii) the code can now be used in the kinetic scheme at arbitrary values of the infrared cutoff 
Ij, and with other definition of the quark masses and of the expectation values. In particular 
the MS definition of rric is now implemented; in) in the kinetic scheme we now normalize the 
matrix element of the Darwin operator, p|,, at the same scale fi as the quark masses and the 
other expectation values (in [S] it was normahzed at /i = 0); iv) the code now computes the 
ratio R*, defined in ^ and necessary to extrapolate the rate measured at the B factories 
to the total semileptonic width. Most of these changes have already been included in the 
version of the Fortran code employed by HFAG |3] in the last few years. 

The 0{al) corrections that are not enhanced by /9o — we will call them non-BLM cor- 
rections — are known to be subdominant when a^ is normalized at m^. They have been 
recently computed in [HlITnillS]- While Refs. [HIITB] adopt numerical methods and can take 
into account arbitrary cuts on the lepton energy, the authors of [I5] expand the moments in 
powers of nic/mb and provide only results without cuts. The two calculations are in good 
agreement and their implementation in our codes is in principle straightforward. However, 
the strong cancellations occurring in the calculation of normalized central moments require 
a high level of numerical precision. Indeed, radiative corrections to the Ei spectrum tend 
to renormalize the tree level spectrum in a nearly constant way, i.e. hard gluon emission is 
comparatively suppressed. This implies that the perturbative corrections tend to drop out 
of normalized moments. Let us consider for instance the first leptonic moment in the kinetic 
scheme with /x = 1 GeV, using r = 0.25, rrih = 4.6 GeV and E^ut = 1 GeV: 



{Ei)E,>iGey = 1.54 GeV 
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It is interesting to note that such kinematic cancellations between numerator and denom- 
inator (identified by the subscript den) affect the 0{as), 0(a^/3o), and two-loop non-BLM 
corrections in a similar way. We have indicated in brackets the numerical uncertainty of 
the non-BLM correction [13]: the resulting coefficient in that case is —0.06 ± 0.12. Similar 
conclusions can be drawn at different values of the cut and for higher linear moments. As 
discussed in [16] , these cancellations are not accidental. In the limit ^ — )■ ^max = 1 — r^ the 
cancellations between numerator and denominator are complete at any perturbative order: 
therefore the higher the cut, the stronger the cancellation. Moreover the peak of the lepton 
energy distribution is relatively narrow and close to the endpoint, which further protects the 
moments from radiative corrections. 

In the case of the higher central moments, additional cancellations occur at each pertur- 
bative order between normalized moments. In £2, for instance, (Ef) and (Ei)'^ tend to cancel 
each other: for the same inputs as in Eq.([6]) we have 

£2 = (Ej) - {Ee)^ = (2.479 - 2.393) GeV^ = 0.087 GeVl 

Such cancellations are quite general and are further enhanced by higher Ecut- They are 



simply a consequence of the fact that, as we have just seen, at each perturbative order the 
spectrum follows approximately the tree-level spectrum, which is peaked at .^ ~ 0.7 — 0.8. 

One obvious consequence of the cancellations we have just discussed is that the numerical 
accuracy with which the non-BLM corrections are known becomes an issue. While the origin 
of the cancellations is understood, we need a precise calculation to know their exact extent, 
and the result will have some impact on the estimate of the remaining theoretical uncertainty. 

The building blocks of the perturbative calculation are the adimensional moments 

Ln = ^— dE,E^— (6) 

where Fq is the total width at the tree-level. Any L„ can be expanded in as and l/rrib 

where we have distinguished between BLM and non-BLM two-loop corrections (/3o = 11 — 
|n;) and indicated by Ln the power-suppressed contributions. The expansion for the 
normalized moments is 
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and an analogous formula holds for R*. 

Both available non-BLM calculations have been performed in the on-shell scheme and give 
US results for L„ at different values of r with an uncertainty due either to the numerical 
integration (and therefore of statistical origin) or to the truncation of the r expansion. 
However, L„ has been computed at ^ 7^ only numerically [16]. The two calculations can 
be combined in order to reduce the final uncertainty. In the analytic calculation of [15] the 
expansion of the 0(0^) corrections to the moments at ^ = includes at most 0{r'^) terms 
and converges quite slowly for r ~ 0.2 — 0.25, in the relevant physical range. We can take 
the size of the last term included, Olr"^), as a rough estimate of its non-gaussian uncertainty. 

It turns then out that the combinations of two-loop non-BLM integrals L) that enter the 

(2) 
normalized moments, 7]l , at ^ = and r<0.26 can be more accurately determined using 

[15] than using the tables of [T6] . 

It is also helpful to notice that, since the electron energy spectrum must vanish at low en- 
ergies at least like Ej, only terms 0{^^) and higher are relevant: the two-loop corrections 7]^ 
must be flat for small ^ and the results of [16] are perfectly consistent with this requirement. 
One can therefore perform a fit to all the available results at different r and ^ values using 
simple functions of ^ and r. We have checked that, given the level of accuracy provided in 
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Figure 1: Combinations of two-loop non-BLM contributions entering the normalized 
leptonic moments and R*: numerical evaluation [16] with blue errors and analytic one 
^ = with red errors vs. fits (shaded bands). 
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[15] at 



J16j . it is sufficient to consider an expansion up to 0{C,^) with coefficients linearly dependent 
on r. The x^/dof is generally low, about 0.4, as the spread of central values in the low-.^ 
region is much smaller than the size of the errors. We also know that the functions 'ql (^, r) 
must vanish linearly at the endpoint ^ = 1 — r^, and one can implement this constraint in 
the fit, but the result are unchanged for ,^ < 0.7. 

In the first plot in Fig. hi we show the combination r][ that enters the first normalized 
leptonic moment ii as a function of ^ for r = 0.25. In the plot we compare the numerical 
evaluations given in Table 1 of [16] and their associated error bars with the fit. At ,^ = we 
also show the result of [IS], whose uncertainty is estimated as explained above. The shaded 
band represents the la uncertainty of the fit. In the case of Eq. (pi) the error of the non-BLM 
0{al) coefficient estimated in this way is ±0.03, with a sizable reduction wrt that equation. 
For higher cuts and for values of r at the edge of the range 0.2 < r < 0.28 considered in 

we show the r-dependence of ^j (^ = 0), comparing 
ear that a linear fit is perfectly adequate. The errors 



16] the reduction is weaker. In Fig. h 
the fit with the existing results. It is c 
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Figure 2: The non-BLM contribution r]\ '{^ = 0) to the first leptonic moment as a function of 
r: numerical evaluation [16J with blue errors and analytic one [15j at ^ = with red errors vs. 
fits (shaded bands). 
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obtained by combining in quadrature the numerical uncertainties of L\ and Lq appear 
overestimated. We have found similar results for the higher linear moments, but we show in 
Figjljonly the combinations that actually enter the second and third central moments, 
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where Wi = L\ /Lq are tree-level functions. Despite the residual uncertainty, the new 
results confirm that the pattern of cancellations observed at 0{as) and 0{alPo) carries on 
at the complete O(a^). This is illustrated in Tables [l] and [2] for Ecut = and IGeV. In these 
Tables we report the values of the first three central moments for the reference values of the 
input parameters 



rUb 



4.6 GeV, 



rric 

pI 



1.15 GeV, 
0.2GeV^ 



pis 



0.4GeV^ 
= -0.15GeV^ 



Pg 



0.35 GeV^ 



(9) 



and for /i = (corresponding to the on-shell scheme) and fi = IGeV. The kinetic scheme 
expressions are obtained from the on-shell ones after re-expressing the pole quark masses 
and their analogues for /x^ and p^ in terms of low energy running quantities, employing the 
two-loop contributions computed in [29] and re-expanding all perturbative series]^ 

In our calculation we remove all terms of 0{asA'^/ml) because they are not yet known 
completely, but we retain suppressed terms of 0{alfi^/ml) that originate in the kinetic 
scheme from a simultaneous perturbative shift in rrif, and /i^ or p^, although they turn out 
to be completely negligible in the case of leptonic moments. 



^Unlike the 0{a^) calculation of the moments, the kinetic scheme expressions of Ref. [29] do not explicitly 
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0.1816(2) 


-0.0340(6) 


tot error [0] 








0.0125 


0.0055 


0.0026 



Table 1: The first three leptonic moments for the reference values of the input parameters and 
-E-cMt = 0, in the on-shell and kinetic schemes. In parentheses the numerical uncertainty of the 
BLM and non-BLM contributions (see text). 
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Table 2: The first three leptonic moments for the reference values of the input parameters and 
Ecut = IGeV, in the on-shell and kinetic schemes. 
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Table 3: The first three leptonic moments for the reference values of the input parameters and 
Ecut = IGeV, in the kinetic scheme with MS charm mass evaluated at /i = 2 and 3GeV, with 
mc(2GeV) = l.lGeV and mc(3GeV) = IGeV. The uncertainty in the 0(«^) is larger in the 
second case because the m^c/m^i, value is closer to the edge of the range considered in [16j. 



Non-BLM contributions to the leptonic moments are generally tiny in the kinetic scheme. 
This is not necessarily the case if we adopt other renormalization schemes. We illustrate the 
point by using an MS definition of the charm quark mass. As mentioned in the Introduction, 
this may be useful for the inclusion in semileptonic fits of precise mass constraints. Table [3] 
shows the results at Ecut = 1 GeV in the case the charm quark mass is renormalized in the 
MS scheme at two different values of the MS scale; jl = 2 and 3 GeV. The bottom mass 
and the OPE parameters are still defined in the kinetic scheme. We observe that in this case 
the non-BLM corrections are not always negligible, especially at /i = 3 GeV. This is due to 
the fact that the MS running is numerically important and driven by non-BLM effects. The 
sizable non-BLM corrections found in this case are therefore related to the change of scheme 
for rric which is known with high accuracy. 

To better understand what drives larger corrections in the MS scheme, let us now consider 
^1,2,3 for Ecut = I GeV and reference inputs. A shift in m,b,c induces in the central leptonic 



include finite charm mass effects in the loops. This mismatch will be ignored in the following: it is effectively 
equivalent to a perturbative redefinition of our mass parameters and it is relevant only when we translate 
them to other schemes, see Section 4. 
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-0.21(30) 


1.2 


-0.30(6) 


-0.32(3) 
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1.5 


-0.33(9) 


-0.31(4) 


-0.27(2) 


+0.3(1.4) 



Table 4: Ratio of non-BLM to BLM contributions to the leptonic moments at various Ecut 
values, in the on-shell scheme. 

monients the shifts 

5£i = 0.325 5mb - 0.237 5mc, 5^2 = 0.094 5mb - 0.058 5mc, 

Sis = 0.009 Srrih - 0.010 5mc, (10) 

which vanish for Snic/Smb ~ 1.4, 1.6 and 0.9, respectively. Therefore, the leptonic central 
moments, and particularly the first two, give similar constraints in the {rric, rrii,) plane. In the 
case of a change of mass scheme, the bulk of the shift in the leptonic moments is given by the 
above 6ii, where Sruc^b represent the difference between the on-shell masses and the masses 
in the new scheme. In the kinetic scheme with /x = 1 GeV we have SrUc/Smb ~ 1.2 at O(a^), 
which explains why the change of scheme does not spoil the cancellations in the perturbative 
corrections, as shown in Table 2. In other schemes and for different values of /i the situation 
can be different, and indeed when we adopt the MS scheme for tUc with ft = 3 GeV (lower 
sector of Table |3]) we observe slightly larger perturbative corrections. As our calculation is 
complete at NNLO, one could worry that higher orders in 6mc,b could spoil the cancellations 
that occur in the on-shell scheme and lead to an underestimate of higher order corrections. 
This appears to be unlikely, except for large MS scales; even that unnatural case, however, 
would be under control, as the 0(0;^) ft evolution is known. 

While in the case of ii and R* the numerical accuracy is always adequate for the fits to 
experimental data, the poor direct knowledge of irreducible non-BLM corrections to the sec- 
ond and especially the third central moment suggests to adopt for them a different approach. 
We have already argued that the two-loop calculations confirm that the cancellations occur- 
ring at 0{as) and 0{a^f3o) replicate at O^a^) as well. This is further illustrated by Table 
[4] where the ratio between the non-BLM and BLM contributions to ?7i,2,3 are shown and, as 
usual, the uncertainty is dominated by that of the non-BLM fits. The ratio of non-BLM to 

(2) 

BLM contributions is well determined and close to —0.3, but in the case of ri2c the uncer- 
tainty is much larger, especially at large cuts, namely where the cancellations are necessarily 
enhanced. We therefore choose to assume a value —0.2 ± 0.1 for this ratio and adopt it 

(2) 

as default. In the case of 773^ the uncertainty is so much bigger than the effect that we 
simply remove this contribution, retaining however reducible non-BLM contributions that 
are important in the MS scheme. Reference values for this default choice are reported in 
Tables [2} [3] in the row denoted by **. 
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Figure 3: /i and Hc dependence of the first three leptonic central moments at Ecut = normalized 
to their reference value fi = fXc = ^ GeV. The three plots refer to £1,2,3, respectively. 



2.1 Scale dependence of leptonic moments 




2.0 2.5 3.0 3.5 

MSbar charm scale (GeV) 



Figure 4: Dependence of the first three leptonic central moments at Ecut = on a) the charm 
MS scale ft (left panel) and b) the scale of a^ when the kinetic scheme is applied to rric as 
well (right panel). The moments are normalized to their values at /2 = 3 GeV and /!„ = rrih, 
respectively. The kinetic cutoff is fixed at 1 GeV. The black, red, blue lines refer to £1,2,3- The 
solid (dashed) lines refer to four (two) loop evolution. 



The size of uncalculated higher order perturbative corrections can be estimated in various 
ways. We have aheady made a few remarks in this direction and we have studied the size 
of the NNLO contributions in various mass schemes. Now we study how our predictions 
for the leptonic moments depend on various unphysical scales that enter the calculation. In 
Fig. |3] we show the dependence of the leptonic moments in the kinetic scheme on the cutoff 
/i and on fic, the cutoff related to the kinetic definition of the charm mass. Indeed, there 
is no reason of principle to set /ic = /i as we have done above. In the plots we have taken 



11 



the fi = fie = ^ GeV values of the b and c masses and of the OPE parameters //^ and p|, in 
Eq. ([9]) and have evolved them to different /x, /^c using 0{a'l) perturbative expansions [29] 
with as{mf,). The plots refer to Ecut = 0, but the absolute size of the scale dependence 
is similar at different cuts, while the relative size varies in an obvious way. Notice that for 
values of fi significantly higher than 1 GeV the kinetic definition may lead to artificially large 
perturbative corrections, especially in the case of the charm mass. In the acceptable range 
shown in the plots the scale dependence is small, and only for ii it is larger than the whole 
O(a^) correction. 

The plot on the left of Fig. |4] refers instead to the calculation with the charm mass 
defined in the MS. It shows the dependence on the charm mass scale /i for fixed kinetic 
cutoff ^ = 1 GeV, normalized to the values fi = 3 GeV, r7ic(3GeV) = 1 GeV. We have 
evolved Trie using the four loop Renormalization Group Evolution (RGE) implemented in 
[30] . It is clear that in this case higher order corrections present in the four loop RGE (and 
to a lesser extent even in the two loop RGE, see Fig. ^ lead to very sizable corrections and 
that scales below 2 GeV should be avoided. Even above 2 GeV the observed scale dependence 
of ii is larger than its whole perturbative contribution. This is not an indication of large 
higher order corrections in the calculation of the moments, but only in the MS evolution of 
rric. The best procedure to extract rric is therefore to compute moments at a scale /i ~ 2 GeV 
that tends to minimize the corrections to the moments, and then to evolve to other scales 
using the most accurate RGE formula. 

Finally, in the right-hand plot of Fig. |4] we show the dependence of the kinetic scheme 
results on the MS scale of a^, fia, for fixed cutoff /i = /ic = 1 GeV. The evolution of a^ is 
again computed at 4 loops, using ^^(m;,) = 0.22 as input. Due to its direct connection with 
the size of the perturbative contributions, the //„ dependence is generally tiny, with the only 
exception of £3. For instance, using as{mb/2) instead of as{mb) in Table 1, we would have 
ii = 1.3870 GeV in place of 1.3865 GeV. 

We will briefly come back to the scheme dependence of the leptonic moments in Sec. 4. 
Before turning to the hadronic moments, however, it is useful to spend a few words on the 
estimate of the total theoretical uncertainty of the leptonic moments, a subject of great 
importance for the global fits from which we extract \Vcb\- In Ref. ^ and in the fits based 
on it, the overall theoretical uncertainty related to ii was computed by adding in quadrature 
various contributions: the shits in ii due to a ±20% (±30%) variation in /i^ and fi^ (p|, and 
pI^), a ±20MeV variation in mc,b, and a ±0.04 variation in as- The uncertainty estimated in 
this way is reported in the last row of Tables [ij and |2} in all cases it is about five times larger 
than the current experimental one. Following the inclusion of the NNLO calculation we 
can now slightly reduce this estimate. Keeping in mind the remaining sources of theoretical 
uncertaintMj, one could use the same method with a ±10 MeV variation in m^fe and ±0.02 
in as] in this way the total theoretical uncertainty in ii is reduced by ~ 25% wrt the last 
row in Tables [l] and [21 while the improvement is only minor for £2,3- 



^While the 0{asfJ,^/ml) contributions to ii are tiny [THj, the 0(l/mj''^) corrections estimated in Ref. [T7] 
are almost as large as the total error in the last row of Tables [Tl and [2| 
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3 Hadronic mass moments at NNLO 

At the time Ref. [6] was published not even all the 0(0;^) corrections to the hadronic moments 
were available for ^ 7^ 0, and therefore the original calculation in the kinetic scheme was 
rather incomplete. The original code underwent a subsequent upgrade to include the full 
0{(3oal), computed numerically following ^0\ and implemented through precise interpolation 
formulas, and in this form has been employed in recent HFAG analyses. Here we discuss the 
implementation of the complete NNLO corrections and present a few results. The code for 
the hadronic moments shares the same features of the leptonic code, including flexibility in 
the choice of the scheme and of all the relevant scales. 

We recall that due to bound state effects the invariant mass Mx of the hadronic system, 
which is measured experimentally, is related to the invariant mass rrix and energy Cx of the 
partonic (quarks and gluons) final state by 

^i = ^'+2e,A + A', (11) 

where A = Mb — rrif,. Since A ^ O.TGeV is relatively large, we do not expand in A/rrif, and 



retain all powers of A. It follows from Eq. ( 11 ) that the Mx moments are linear combinations 
of the moments of m^. and Cx, namely of the building blocks 

M^j = ^^ / {ml - mlY ei dT . (12) 

For instance, in the case of the first normalized moment we have 

o V o -1-2 o iVfin -;- Mm 

Moo ^mo 

Like in the case of leptonic moments, we consider only central higher moments, and specifi- 
cally 

h, = {MJ,), h, = {{Ml-{Ml)f), h, = {{Ml - {Ml)f) . (13) 

The calculation of /ii,2,3 requires the knowledge of the building blocks M^- with i + j < 
1,2,3, respectively. They receive both perturbative and non-perturbative power-suppressed 
contributions: 



as{mb) - .(1) , fa. 



M,, = Mf + ^^M« + (^) (/3oM,r^) + M^^) + M^"'^ + ... (14) 

where the tree-level contributions M^ • vanish for i > and we have distinguished between 
BLM and non-BLM two-loop corrections (/3o = 11 — |n/, rii = 3). Since we are interested in 
the normalized moments, we will be mostly concerned with the combinations 
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Second hadr energy moment 




0.2 0.3 0.4 0.5 
^ {mjmt=0.22) 
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Figure 5: Combinations of two-loop non-BLM contributions Xoi < X02 . 2nd Xio entering the 
hadronic moments: numerical evaluation [16J with blue errors and analytic one ^5j at ^ = with 
red errors vs. fits (shaded bands). Xio (lower panel) is shown for r = 0.2,0.22,0.24,0.26,0.28. 



with a = 1, BLM, 2, pow, which enter the perturbative and non-perturbative expansion for 
the normahzed moments and are functions of ^ and r: 

M^? _ „ _„ 

(16) 



M,, 



M 



00 



M, 



(0) 

00 



~r Xij ~r 



TT 



/3ox!"^"^ 



1,2 J 



+ Xij 



lAij 



M, 



(0)- 

00 



+x!r^+... 



The first 0{as) calculation appeared in [31] and was later completed in [321 |33l HD] , while 
the 0{a1f3o) (BLM) calculation was presented in [331 HD]. Analytic expressions for the 0{as) 
contributions to the building blocks Mij with z 7^ are available as functions of ^,r [TO] . 
For the remaining building blocks and for all the 0(a^/3o) we employ high-precision two- 
parameters interpolation formulas to the results of the numerical code of Ref. [TU], which 
are valid in the range < .^ < 0.8 and 0.18 < r < 0.29. The use of these approximations 
induces a negligible error. 

The complete NNLO results presented here are based on the on-shell scheme results of 
Refs. [151 CH]- R-sf. [15] provides us with NNLO contributions to the first two e^ moments 
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E. 



cut 



PI 

Xoi 
(2,BLM) 
Xoi 



[2) 
Xi(/ 
(2,BLM) 
XlQ 



(21 
Xo2 
(2.BLM) 
X(12' 



PI 
Xll 
(2.BLM) 
Xll 



(2) 
Xi2 
(2, ELM) 
Xi2 





0.6 
0.9 

1.2 
1.5 



-0.24(1) -0.24 -0.23(1) -0.24 -0.24 

-0.25(2) -0.24 -0.24(1) -0.24 -0.24 

-0.26(3) -0.23 -0.24(1) -0.23 -0.24 

-0.26(5) -0.23 -0.24(2) -0.23 -0.23 

-0.14(20) -0.23 -0.17(7) -0.23 -0.24 



Table 5: Ratio of various perturbative contributions relevant to the hadronic moments at various 
Ecut values, in the on-shell scheme. 





/i = 


H = IGeV 




hi 


h2 


hs 


hi 


h2 


hs 


LO 


4.420 


0.199 


-0.03 


4.420 


0.199 


-0.03 


power 


4.515 


0.767 


5.51 


4.515 


0.767 


5.51 


0(a.) 


4.662 


1.239 


7.78 


4.504 


1.069 


6.13 


0(/3oaD 


4.834 


1.661 


9.51 


4.493 


1.301 


6.39 


0{a's) 


4.811 


1.603(5) 


9.22(7) 


4.495 


1.222(5) 


6.24(7) 


tot error |6] 








0.143 


0.510 


1.38 



r(2) nA2) 



Table 6: The first three hadronic moments for the reference values of the input parameters and 
Ecut = 0, in the on-shell and kinetic schemes. 

at ^ = 0, and we have the non-BLM contributions to several of the building blocks from 
Ref. [16]. In all the available cases we proceed in the same way as for leptonic moments and 
perform fits to the tables of (TB] and to the ^ = of [T3]; some of the results are shown 
in Fig. p] Unfortunately not all the building blocks needed for the second and third M 
moments are reported in the tables of [I6]: the missing building blocks are M; 

(2) 

MgQ , which have been computed only for r = 0.25 and for .^ = and 0.435 

In order to deal with these three cases we observe that in the on-shell scheme the ratio 
of non-BLM to ELM contributions to Xij is always remarkably insensitive to the value of 
^, as shown in Table [5| and in fact much more insensitive than the analogous ratios to the 
tree-level or one-loop contributions. In particular, this ratio is known precisely for z 7^ and 
varies by about ±0.01 in the range 0.2 < r < 0.28. In view of the other sources of theoretical 
errors, it is therefore sufficient to estimate the missing X20 > X30 ; X21 using the values of their 
ratios in the two available points, and to assign them a ±0.02 error: 



20 ) ^21') and 



(2) 
X20 
(BLM) 
A 20 



-0.215±0.020, 



(2) 
X2I 
(BLM) 
A21 



-0.215±0.020, 



(2) 
X30 
(BLM) 
A 30 



-0.205±0.020 (17) 



As illustrated in Fig. [5] by a few examples, the precision of the fits to the functions 
Xoi is similar to that of the functions 77) entering the leptonic moments. On the other 
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/i = 


/i = IGeV 




hi 


h2 


hz 


hi 


h2 


h. 


LO 


4.345 


0.198 


-0.02 


4.345 


0.198 


-0.02 


l/ml 


4.452 


0.515 


4.90 


4.452 


0.515 


4.90 


0{as) 


4.563 


0.814 


5.96 


4.426 


0.723 


4.50 


0{P,al) 


4.701 


1.105 


6.85 


4.404 


0.894 


4.08 


0{al) 


4.682(1) 


1.066(3) 


6.69(4) 


4.411(1) 


0.832(4) 


4.08(4) 


tot error ^Gj 








0.149 


0.501 


1.20 



Table 7: The first three hadronic moments for the reference values of the input parameters and 
Ecut = IGeV, in the on-shell and kinetic schemes. 





/i = IGeV, m^s^2GeV) 


/i = IGeV, mf^i'iGeV) 




hi /l2 /i3 


hi /l2 /l3 


l/ml 


4.301 0.551 4.94 


4.020 0.618 5.02 


0(a.) 


4.355 0.758 4.60 


4.192 0.830 4.79 


0{^al) 


4.304 0.936 4.21 


4.169 1.015 4.49 


0{al) 


4.328 0.865(4) 4.18(4) 


4.245(1) 0.922(5) 4.38(4) 



Table 8: The first three hadronic moments for the reference values of the input parameters and 
Ecut = IGeV, in the kinetic scheme with MS charm mass evaluated at /i = 2 and SGeV, with 
mc(2GeV) = l.lGeV and mc(3GeV) = IGeV. The uncertainty in the 0(0"^) is larger in the 
second case because the m^c/m.^, value is closer to the edge of the range considered in [16J . 



hand, the numerical evaluation of \i6l is very accurate in the case of xlj for ^ > and a 
quadratic dependence on r has to be included in the functional form. As the cancellations 
between different contributions are not as severe as for leptonic moments, the accuracy of 
the non-BLM corrections is in fact a minor issue in the case of hadronic moments. 

Numerically, the non-BLM corrections in the on-shell scheme tend to partly compensate 
the ELM corrections and, as perturbative corrections in general, are more important than in 
the leptonic case, see Tables |6] and [TJ The non-BLM effect is about 0.4%, 4%, 3% for /ii,2,3, 
respectively. In the kinetic scheme with /i = 1 GeV the non-BLM corrections to the first and 
third moment are slightly suppressed, while the non-BLM correction to /i2 is larger than in 
the on-shell scheme (7-9%), see also Fig. [6] It is worth reminding that the perturbative con- 
tributions to the building blocks Mu, M20 and Mio (all vanishing at tree-level) are dominant 
in /i2; and that these building blocks are known only at the next-to- leading order level. This 
explains the large relative value of the non-BLM corrections and implies a commensurate 
uncertainty. But the results shown in the Tables hide an important cancellation between the 
/i^ and p|) contributions to /i2 for the reference values of Eq. (|9]): for Ecut = 1 GeV they are 
about +1.6 and — 1.2GeV^, respectively. Non-perturbative contributions indeed dominate 
the higher moments /12 3, which are therefore subject to much larger uncertainties from both 
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Figure 6: E^ut dependence of leptonic (left) and hadronic (right) at NNLO in the kinetic scheme 
with /i = /ic = 1 GeV. The black, red, blue lines refer to £i, 10x^2, —10x^3 and to /ii, 8x /i2, /^s, 
respectively, each expressed in GeV to the appropriate power. The dotted lines (indistinguishable 
for the leptonic moments) represent the predictions at 0(a^/3o). 



higher order terms in the OPE and the missing 0[ol^ corrections to the Wilson coefficients 
of the dimension 5 and 6 operators. 

3.1 Scale dependence of hadronic moments 

The [I and /i^ dependence of the hadronic moments in the kinetic scheme is shown in Fig. [7] 
for Ecut = 1 GeV. The first moment has very small residual scale dependence, like the 
leptonic moments. Since the second and third hadronic moments receive much larger power 
and perturbative contributions, the larger /x and especially /Xc dependence is not surprising. 
The results of the calculation of the moments when the charm mass is renormalized in 
the MS scheme are illustrated in Table [SJ for two different values of the charm scale Jx. It 
is interesting that for /i = 3 GeV the change of scheme induces large non-BLM corrections, 
suggesting that /i = 2 GeV might be a better choice. Fig. |8] is the analogue of Fig. |4] for 
the hadronic moments. The plot on the left shows the dependence on the charm mass scale 
/2, when /ii,2,3 are normalized to their values for fl = 3 GeV, mc(3GeV) = 1 GeV. In this 
case the /2 dependence is even larger than in the leptonic case, but again it should not be 
interpreted as indication of large higher order corrections in the calculation of the moments 
and one should preferably use values fl that minimize the corrections to the moments. 

Like the leptonic moments, the first hadronic moment can be used to constrain the 
heavy quark masses. The higher hadronic moments are too sensitive to higher dimensional 
expectation values to play a role in this respect. At reference values of the inputs and for 
Ecut = 1 GeV, a small shift in rrib c induces 



5hi 



-4.84 5mb + 3.08 (5m, 



which vanish for 6mc/Smb ~ 1.6. Therefore, hi gives a constraint in the {mc,mb) plane very 
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Figure 7: jj, and fic dependence of the first three hadronic central moments at E^ut = 1 GeV 
normalized to their reference value fx = f^c = ^ GeV. The three plots refer to /ii,2,3, respectively. 



similar to the leptonic central moments. The considerations on scheme dependence made 



after Eq. (10) apply to the first hadronic moment as well. 



Finally, in the right-hand plot of Fig. [8] we show the dependence of the kinetic scheme 
results on the MS scale of as, fia, for fixed cutoff /x = /ic = 1 GeV. The evolution of as is 
again computed at 4 loops, using as{mb) = 0.22 as input. Due to its direct connection with 
the size of the perturbative contributions, the //„ dependence is small for hi. The strong 
residual dependence of /12, on the other hand, suggests that in the absence of higher order 
corrections, it would be prudent to choose /Xq close to mb/2. 

The overall theoretical errors computed using the method proposed in p] are reported 
in the last row of Tables [6] and [7j they are generally a factor two or more larger than the 
present experimental ones. Cutting by half the variation in rrib^c and a^, as it was proposed 
in the previous Section, leads to a 25-30% smaller total uncertainty on hi, while the error on 
/i2,3 is almost unaffected. The calculation of 0{as) corrections to the Wilson coefficients and 
that of higher orders in the OPE |17|, will have a stronger impact on the total uncertainty. 
We also recall that the 0(l/r7i^' ) contributions to hi found in Ref. [l7j are of the same order 



of magnitude of the total error in the last row of Tables^ and ^ while those of 0{asiJ''i) [18 
are comparatively small. 



4 Mass scheme conversion 

Semileptonic moments provide interesting constraints on the bottom and charm masses, 
which can be compared and combined with other experimental determinations of these pa- 
rameters. Since in many applications the quark masses are renormalized in the MS, we will 
now review the relation between the kinetic and MS definitions of the heavy quark masses. 
Accurate conversion formulas can also be used to study the scheme dependence of our re- 
sults. The perturbative relation between the quark pole mass and the kinetic mass is known 
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Figure 8: Dependence of /ii,2,3 at Ecut = on a) the charm MS scale ft (left panel) and b) the 
scale of as when the kinetic scheme is applied to rric as well (right panel). The moments are 
normalized to their values at /i = 3 GeV and jja = f^b, respectively. The kinetic cutoff is fixed 
at 1 GeV. The black, red, blue lines refer to /ii,2,3- 



completely at 0{a1) and 0(q;^/9o) IZS]- Combining it with the one between pole and MS 
masses, known to 0{al) [33], one gets the following NNLO expansion for the ratio between 
the MS mass fn{fl) and its kinetic counterpart m{fi) 
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where we have used 
L = In [/iVm(^)2] 



/3t^ 



11 - 2/3 n;, A(x) 



TT 



-X - 0.597x^ + 0.230a;^ 



The function A{x) arises from the two-loop diagrams with a closed massive quark loop in 
the relation between pole and MS mass [351 EO]- It represents the correction to the massless 
loop and is only relevant in the case of charm mass effects on the b mass. As mentioned in 
footnote [ij such effects are not included in the conversion to the kinetic scheme, which is 
performed assuming charm to be decoupled |29Fj There is also a dependence on mh, the 



■^Finite charm mass 0{a^) effects in the transition between pole and kinetic mass can in principle com- 
puted from Eq. (11) of ^. 
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mass scale that defines the threshold between rii and rii + l quark flavors, generally set equal 
to m. In the limit /j, ^ the above formula relates the pole and the MS masses. 

When the scale of the MS mass is equal to the MS mass itself and the coupling constant 
is evaluated at the mass scale the above formula reduces to 
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(19) 



This result differs from those reported in [221 E] by a term (^) ( is ~'~ ^ ) ■ '^^^ difference 

can be traced back to a mismatch between the number of quarks considered in the MS-pole 
relation. In the limit /x — )• 0, Eq. (19) reproduces the standard two-loop MS-pole relation, 

have been also considered in 



cf. Eq. 16 of [30] and [36]. The BLM corrections of 0{a^J^ 



[221 E]. 



4.1 Charm mass conversion 

In the case of the charm mass we employ the two above formulas with ni = 3, neglecting all 



A functions. Eq. (19) with a 



{rric) = 0.41 and mc(l GeV) = 1.15 gives 



m. 



(m^) = 1.150 + 0.108^^+ 0.071^2 = 1.329 GeV ^ 



(20) 



where the one and two-loop contributions are identified by the subscript. Unsurprisingly, the 
perturbative series converges poorly. One easily verifies that BLM corrections dominate the 
two- loop contribution and that the result depends strongly on the scale of as- The apparent 
convergence of the perturbative series can be drastically improved if one normalizes the charm 
mass at higher MS scales. In fact, in the evolution of m^(/2) to low scales the effect of RGE 
resummation becomes significant, but these contributions are absent from the fixed order 
NNLO calculation of the moments from which the kinetic mass could be extracted. In order 
to convert the results of a fit to the semileptonic moments in the kinetic scheme into the MS 
scheme (or vice versa) it is therefore important to adopt a scale fi above 2 GeV. To illustrate 
the point let us use ft = 3 GeV and adopt the central value of [23], m^(3GeV) = 0.986 GeV. 



„(3) 



Inverting Eq. (18) and employing al (3GeV) = 0.247 (which corresponds to al {rrib) = 0.22 



(4)/ 



for irih = 1.5 GeV), we obtain the kinetic charm mass 



mf"(lGeV) 



0.986 + 0.083„, + 0.022, 



1.091 GeV 



(21) 
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where the final result depends little on the scale of a^ and on its exact value. The difference 

me(l GeV) - m^(3GeV) = (0.105 ± 0.015) GeV (22) 

depends of course on the input m^(3GeV): for values between 0.95 and 1.05 GeV the central 
value moves from 0.096 to 0.121 GeV. 

Using a lower scale /i = 2 GeV the difference between kinetic and MS is even smaller. 
Using m;:(2GeV) = 1.092 GeV, obtained from the four-loop evolution of the result of jS3] . 
one finds 

m^*"(lGeV) = 1.092 + 0.037„, - 0.013^2 = 1.116 GeV (23) 

with very small dependence on the scale of as- The difference between the final results in 



Eqs. (21 ) and (23) is mostly due to our use of four-loop RGE to connect the inputs at /i = 3 



and 2 GeV. Our estimate is 

m,(l GeV) - m^(2GeV) = (0.024 ± 0.010) GeV, (24) 

whose central value varies between 0.017 and 0.032 GeV for m^(2GeV) in the range 1.05 to 
1.15 GeV. Eqs. ( 22p4 ) can be used to compare accurately the predictions for the semilep- 



tonic moments when the kinetic and MS schemes are used for mc. It turns out that the 
size of the scheme dependence found in this way is similar to the uncertainty due to scale 
dependence investigated in Sees. 2 and 3. 

Once the relation between kinetic and MS masses is known at large /2 one can use the four- 
loop RGE together with Eq. (21p3) to compute fn^ifn^). Using the more precise Eq. (23) 



we obtain 

mc(l GeV) - m-c^irrTc) = 0.16 ± 0.02 GeV. 

4.2 Bottom mass conversion 

In the case of the relation between MS and pole bottom mass the effects related to the 
charm mass in closed quark loops and described by A(r) are small but not negligible. For 
realistic values of the charm to bottom mass ratio A(r) compensates up to 40% of the light 



quark contribution, i.e. the term multiplied by rii in the fourth line of Eq. (18). Since the 
corresponding effects are not included in the definition of kinetic m;,, in applying Eqs. ( T8p9 ) 
to the bottom mass we therefore have the following options 

(a) we decouple charm completely and set rii = 3, A(r) = 0; 

(6) we treat charm as massless and set rii = 4, A(r) = 0; 

(c) we decouple charm in the kinetic part of the calculation only, set rii = 3, and replace 



A(r) -> A(r) 



jn_ 

24 192 ■ 



It is clear that the correct result lies somewhere between (a) and (6). On the other hand, since 
the NNLO expressions of the semileptonic moments have been derived under the assumption 
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Figure 9: M dependence of mh(mi,) in case (c) for mi,{l GeV) 
curves refer to one-loop, BLM, two-loop, respectively. 



4.55 GeV. The blue, red, black 





a,{mb) 


mb(l GeV) 


0.21 0.22 0.23 


4.500 
4.550 
4.600 


4.150 4.127 4.103 
4.194 4.171 4.147 
4.239 4.216 4.192 



Table 9: mhlmh) for different values of the kinetic mass and of as{mi,) in case (c) with 
M = 2.5 GeV. 

of charm decoupling in the relation between pole and kinetic masses, but we have otherwise 
kept the charm mass effects, option (c) is the most appropriate to convert the kinetic mass 
extracted from a NNLO fit to semileptonic moments into the MS scheme. 

Let us now illustrate the three options computing the MS bottom mass from the kinetic 
one at a reference value /i = 1 GeV. We use mfe(l GeV) = 4.55 GeV and as{mh) = 0.22 



mirrTb) = 4.550 - 0.290^,, - 0.073^2 
4.550 - 0.290„, - 0.060^2 
4.550 - 0.290„. - 0.057„2 



4. 187 GeV, 
4.200 GeV, 
4.202 GeV, 



(b) 



(25) 



from which we see that option (c) is numerically very close to (b). The BLM corrections 
are large in Eq. (25) but partly (50%) compensated by sizable non-BLM ones. The 0{alP^) 
corrections [29| shift rribirfTi) further down by about 20 MeV. The dependence of the NNLO 
result on the scale of a^ is shown in Fig. [9] it suggests a residual uncertainty of about 40 
MeV. Table [9] shows the value of fnbirfTb) for different values of the inputs, using for a^ a 
lower scale than nib itself, M = 2.5 GeV. The dependence on r is negligible. 

If one computes rrib from m{,(l GeV) at a scale around 3 GeV the perturbative corrections 
are smaller because the numerical values of the masses get closer. For instance one gets 
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m^(2.8 GeV) = 4.550 + 0.023^^ - 0.024„2 = 4.550. This result can be evolved up using 
four-loop RGE, finding fnh{fnb) = 4.182, which is consistent with the values in Table |9] but 
has arguably a slightly smaller error. Our final estimate is 

mb(l GeV) - fniim-b) = (0.37 ± 0.03) GeV, 

The error also includes in quadrature the parametric as error, Sasirrib) = 0.007. 

5 Summary 

I have presented the results of a complete O(a^) calculation of the moments of the lepton 
energy and hadronic invariant mass distributions in semileptonic B — )■ Xciiy decays. The 
numerical Fortran codq^ that has been developed allows us to work in any perturbative 
scheme; I have adopted the kinetic scheme, but have also considered the option of the MS 
scheme for the definition of the charm mass. The non-BLM corrections are generally small 
and within the expected range. The new code will allow for a NNLO fit to the semileptonic 
moments, with the possible inclusion of precise mass constraints, in a variety of perturbative 
schemes. To facilitate the inclusion of mass constraints obtained in the MS scheme, I have 
critically reexamined the conversion from that scheme to the kinetic scheme for the heavy 
quark masses. Additional details concerning the NNLO prediction of the total rate can be 
found in the Appendix. 

I have performed a detailed study of the dependence of the results on various unphysical 
scales, obtaining useful indications on the size of the residual perturbative uncertainty in the 
moments, which is now definitely smaller than the uncertainty due to higher dimensional 
contributions in the OPE and to 0{as) contributions to the Wilson coefficients of dimension 
5 and 6 operators. Fortunately, work on both these aspects is progressing and there are good 
prospects for a more accurate prediction of the inclusive semileptonic moments and a more 
precise determination of \Vcb\- 
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Appendix 

The determination of \Vcb\ from inclusive semileptonic B decays is based on the expression 
for the charmed total semileptonic width. In this Appendix we update the numerical analysis 
in the kinetic scheme given in Ref. [S] with the inclusion of the exact two-loop correction 
computed in [TH [15], HE] . Keeping only terms through 0(1/ ml), the general structure of the 
width is (see [5] and refs. therein) 



r|B^.Y,eP] = ^^K,\'g{p){l + A, 






^^l +(l_ 2(l-pr ^ ^'G-'^ ^ d{p) Pi 



(26) 



32p3 
3 



2ml V2 g{p) J ml g{p)ml 

where p = (mjmb)^, g{p) = 1 - 8p + 8p^ - p"^ - Up'^lnp, d{p) = 81np - i|^ + 

8p^ 3^ + y, and all the masses and OPE parameters are defined in the kinetic scheme 

with cutoff p. The expression for the width is /i-independent through 0{a1). The term 
Af.^ = 2a/TT In Mz/m,h ~ 0.014 is due to the electromagnetic running of the four-fermion 
operator from the weak to the b scale and represents the leading electroweak correction [37j . 
The results of the two-loop BLM and non-BLM perturbative corrections [101 [El [16] in the 
on-shell scheme are well approximated in the relevant mass range by a simple interpolation 
formula, valid for 0.17 < r < 0.3, 

p(2)(r,0) = (-3.381 + 7.15r-5.18r2)/3o + 7.51 -21.46r + 19.8r2. 

Applying these results and the well-known one-loop contribution [39] to the kinetic scheme 
calculation with p = 1 GeV, mb = 4.6 GeV, and r = 0.25, the NNLO perturbative series 
becomes 

r[5^Xcez/] ocl-0.96 — -0.48/3o( — I +0.811 — ) ^0.916 (27) 

where /3o = 9 corresponds to three light flavors, and we have used ^^(mfe) = 0.22 for the nu- 
merical evaluation. For comparison, Ref. [5] used the same inputs and had —1.0 as coefficient 



of al and the global perturbative factor of Eq. (|27j) was 0.908. Of course, the difference is 
due to the recent complete NNLO calculation of [El [15]. As already discussed, finite charm 
quark mass effects in the NNLO conversion to the kinetic scheme have not been included. 
They can be expected to decrease the above value by up to 0.002. Higher order BLM correc- 
tions are also known ^38j and have been studied in the kinetic scheme where the resummed 
BLM result is numerically very close to the NNLO one [5]. The residual dependence of the 
overall perturbative factor on the scale of as is mild: changing it between m,^ and m,b/2 de- 
creases by less than 1%, which represents a reasonable estimate for the residual perturbative 
uncertainty. 
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We also provide an approximate formula for \Vcb\ in terms of the OPE parameters and 
as- Using tb = 1.582 ps and BR^^^c = 0.105 (the measured semileptonic branching ratio 
after subtraction of 6 — ?■ ulu transitions) we have 

-1^ = l + 0.275(a, -0.22) -0.65 (m^- 4.6) + 0.39 (m,- 1.15) + 0.014 (/i^- 0.4) 
0.0416 

+0.055 {fil - 0.35) + 0.10 (p^ - 0.2) - 0.012 [pig + 0.15), (28) 

where all dimensionful quantities are expressed in GeV to the appropriate power and are 
understood in the kinetic scheme with /i = 1 GeV. The numerical accuracy of this formula 
is limited to 1% in the la range of the current HFAG fits [5]. The combination of rricfi in 
the above expression is very similar to the one that appears in the leptonic moments, see 



Eq. (10). 



If the MS scheme is adopted for the charm mass with mc{2GeV) = 1.12 GeV and 
mfe(l GeV) = 4.6 GeV, the NNLO perturbative series expressed in terms of asimi,) becomes 

r[fi^X,ez/] ocl-1.27— -0.33/3o ( — ) -0.29 ( — ) f^ 0.895. (29) 

A simple test of the residual scheme dependence consists in comparing the products of the 
function g{p) with the overall perturbative factor, Eqs. (27p9), computed in different schemes 



for rric- This product should be scheme independent when the same bottom mass is employed. 
In the case considered here the values mc(l GeV) = 1.15 GeV and mc(2GeV) = 1.12 GeV 



satisfy Eq.(24) and the product is scheme independent to excellent approximation. For a 
larger MS scale /2 = 3 GeV the perturbative expansion converges slower, but the scheme 
dependence is well within the 1% uncertainty estimate. 
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